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ABSTRACT 

The spectra of many X-ray pulsars show, in addition to a power law, a low-energy component that 
has often been modeled as a blackbody with fcTeB ~ 0.1 keV. However the physical origin of this soft 
excess has remained a mystery. We examine a sample of well-studied, bright X-ray pulsars, which 
have been observed using ROSAT, ASCA, Ginga, RXTE, BeppoSAX, Chandra, and XMM-Newton. 
In particular, we consider the Magellanic Cloud pulsars SMC X-l, LMC X-4, XTE J0111. 2-7317, and 
RX J0059.2-7138 and the Galactic sources Her X-l, 4U 1626-67, Cen X-3, and Vela X-l. We show 
that the soft excess is a very common if not ubiquitous feature intrinsic to X-ray pulsars. We evaluate 
several possible mechanisms for the soft emission, using theoretical arguments as well as observational 
clues such as spectral shapes, eclipses, pulsations of the soft component, and superorbital modulation 
of the source flux. We find that reprocessing of hard X-rays from the neutron star by the inner 
region of the accretion disk is the only process that can explain the soft excess in all the pulsars with 
Lx > 10 38 ergs s _1 . Other mechanisms, such as emission from diffuse gas in the system, are important 
in less luminous objects. 

Subject headings: pulsars: general, X-rays: binaries, stars: neutron, accretion, accretion disks 



1. INTRODUCTION 

X-ray binary pulsars (XBPs) have provided a fruitful 
laboratory for high-energy astroph ysics since th eir dis- 
covery in the 1970's (for a review see lNagasell989|) . Many 
XBPs have been observed extensively, and are useful for 
studying a variety of astrophysical processes, including 
accretion onto compact objects, binary interactions, and 
the effects of strong magnetic fields. 

In recent years, studies of XBPs with orbiting obser- 
vatories such as ROSAT, ASCA, Gmga, RXTE, Bep- 
poSAX, Chandra, and XMM-Newton have allowed for 
closer scrutiny of the relevant physics. In particular, 
most bright XBPs outside the Galaxy (for which soft 
X-rays are not heavily absorbed by Galactic gas) have 
shown a marked soft excess in t heir X-ray spe ctra above 
the standard power-law model l)Nagasel 12002'). Authors 
have proposed a variety of simple models to describe 
this soft excess, but there has been considerable debate 
over the origin of this feature in a number of sources. 
For example, the soft component in LMC X-4 has been 
modeled as a the rmal bremsstrahlung-shaped fan beam 
l)Woo et al.ll996j), power-law em ission from the upper ac- 



cretion column (|Paul et all2002fl , diffuse bremsstrahlung 
emission, and/or re processing by accreting material at 
the magnetosphere (|La Barbera et afll2001[h In this pa- 
per we seek to narrow these candidate mechanisms for 
the soft excess, by examining the physical processes in 
light of a wide range of observations. 

The Galactic XBP Hercules X-l has a well-studied soft 
excess that can be explained by reprocessin g of hard X- 
rays by the inner edge of t he accretion disk l|Endo et al.l 
120001 iRamsav et al.ll2002j) . At the moment Her X-l is 
the XBP for which the soft excess is best observed and 
understood, and we will use it to help guide our study of 



other sources. 

Because the Galaxy's disk contains X-ray-absorbing 
gas, the soft excess is not easy to observe in most Galac- 
tic XBPs. This emission is best seen in the Magellenic 
Clouds, at high Galactic latitude where there is not much 
low-energy absorption. The bright XBPs in these nearby 
galaxies are generally more luminous than those in the 
Milky Way and are close enough to obtain good spec- 
tra with a high signal-to-noise ratio. In a survey o f 
ASCA observations of the SMC, lYokogawa et alJ (2003) 
found 30 XBPs, of which three have X-ray luminosities 
large enough to extract good spectra. All three of these 
sources (SMC X-l, RX J0059. 2-7138, and XTE J0111.2- 
7317) show a clear soft excess feature in the spectrum. 
The one comparably bright XBP in the LMC (LMC X-4) 
also shows a soft excess. All four of these sources have 
Lx ~ 10 38 ergs s _1 . 

4U 1626-67, Centaurus X-3, and Vela X-l in the 
Galaxy, and EXO 053109-6609.2 a nd A 0538-66 in the 



LMC also show soft e xcess features ( Angelinij^dJ 
[ Burderi et all 1200(1 iHaberl I fl99l lHaberl et all 
iMavro matakis & Haberll993|) . In addition, a recent sur- 
vey of the SMC with XMM has revea led soft component s 
in several lower luminosity pulsars l|Sasaki et alj f2003). 
In § 2 we summarize observational studies of the XBPs 
mentioned above, and in § 3 we examine how common 
the soft excess feature is in the total population of XBPs 
and introduce possible mechanisms for its origin. In §§ 
4-7 we evaluate these possibilities in terms of both the- 
ory and observations. In § 8 we discuss the implications 
of the results and future work, and in § 9 we summarize 
our results. 
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TABLE 1 

Orbital and spectral parameters for XBPs with soft excesses 
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-•super orbital (d) 

L x (ergs s 1 ) 

L ao ft/Lx 
Soft pulses? 

Nn (10 21 cm" 2 ) 

-^hard 

E C ut (keV) 
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2. OBSERVATIONAL BACKGROUND 

Here we review a number of observations of XBPs with 
a soft excess. Orbital and spectral parameters for each 
source are compiled in Table 1. 

2.1. Her X-l 

Her X-l was one of the first known X-ray pulsars 
{Tananbaum et aTlll972|) and is one of the few with a 
low-mass companion. It has a pulse period of 1.24 s 
and an orbital period of 1.7 days. The system is at 
high inclination, and the A7 companion star, with a 
mass of ~2 Mp,, eclipses the pulsar every orbital period 
{Middleditch fc Nelsonll 1976ft . In addition there is a 35- 
day cycle in the X-ray intensity, consisting of "main-on" 
and "short-on" states separated by two "low" states in 
which the X-ray flux drops b y several orders of magni- 
tude { Jones &: Formanl Il976|) . This modulation is be- 
lieved to be due to a tilted precessing accretion disk 
that periodically blo cks the central X-ray source (see 
iScott fc Leahvlll999l for a review). Given a distance of 
5 kpc, Her X-l has L x ~2x 10 3 '' ergs s -1 . 

The spectrum of Her X-l has been extensively studied, 
and is generally fitted with a cutoff power law with T ~ 1, 
a blackbody with /cTbb ~ 0.1 keV, and often broad-line 
emission at ~1 keV. There is also iron emission at ^6.5 
keV a nd a cyclotron feature at ~40 keV {dal Fiume et al.l 
1998). The energy of the (fundamental) cyclotron line is 
related to the B field by 

E cyc = 11.6(5/10 12 G)(l + z)- 1 keV, 

wher e z is the gravitational redshift fi.e.. ICoburn et ahl 
2002). For surface emission from a neutron star with 
R ns = 10 km and M = 1.4M©, we have z ~ 
0.3 {Coburn et alJ I2002T) . Usually neither E cyc nor z 
is well constrained, so for simplicity we take B ~ 
10 12 (£ cyc /10 keV) G in this paper. Thus, Her X-l has 
B~4x 10 12 G. 

The pulsations of the soft component are generally 
smoo ther than and out of phase with the hard compo- 
nent l|McCrav et al.lll982tlEndo et alJl2000h. This phase 
offset varies over the superorbital period jllamsav et al.l 
2002) , supporting the idea that the soft excess is emitted 
by the precessing accretion disk. We discuss this process 
in detail in § 7. 

2.2. The luminous Magellanic Cloud XBPs 

SMC X-l is a luminous, eclipsing, high-mass XBP sys- 
tem. It was fir st identified as a po int source in Uhuru 
observatio ns by |Leong e t alJ i!971|) . Eclipses were dis- 
co vered bv | gchreier et alJ {1972al) and X-ray pulsations 
by iLucke et al.1 {19761 ) The optical count erpart is Sk 
160, a B0 I supergiant {Webster et al Jl972t lLilieillT97^ . 
The orbit is close to circular, with orbital period ~3.9 
days and pulse period 0.7 s. The source has a high 
Lx ~ 2 x 10 38 ergs s _1 . As in Her X-l there are large, 
quasi-periodic variations of an ~60 day period between 
high and low states of the X -ray flux, attributed to block - 
ing by a precessing disk {Oruber fc Rothschild! (19841. 
These variations are stable, and the period of the cy- 
cles o scillates from ~50 to ~60 days {Woidowski et all 
1998). This oscillation is it self periodic on the order of 
years {Clarkson et alJ l2003). The spectrum of SMC X-l 
is usually fitted as a cutoff power law with T ~ 0.9, an 



iron line at ^6.4 keV, and a soft component modeled as 
a broken power law, therma l bremsstrahlung, or a black- 
body with kT BB ~ 0.16 keV. lPaul etHTI l|2002j) fitted four 
different models to the soft component and found that all 
gave satisfactory fits. Similarly to Her X-l, they found 
that the soft component pulses roughly sinusoidally, out 
of phase with the sharper-peaked power law. 

LMC X-4 is a high-mass XBP very similar to SMC 
X-l, and was also d iscovered by Uhuru observations 
(iGiaccon i et al.l I 1972D. Its co mpanion is an III-IV 
star (iSanduleak fc Philip|ll977|) . It is eclipsing {Li et alJ 
119781 iWhitel 119781) and has a pulse period of 13.5 s 
(Kclle vet al.lll983|) and an orbital period of 1.4 days. 
LMC X-4 also has a stable superorbital pe riod of ~30 
days, again attributed to a precessing disk {Lang et alJ 
1981). About once per day the source shows flares lasting 
^20-45 minutes during which the X-ray fl ux increases by 
up to -20 times fe.g.. lLevine et al.lll991|) . Like SMC X- 
1, the high-state spectrum is fitted by a cutoff power law 
(F ~ 0.6), with a soft component modeled in a variety 
of ways suggesting different emission mecha nisms. In the 
low s tate, the power law hardens to T ~ 0.1 {Naik fe Paull 
2002). There are iron emis sion at ~ 6.5 keV and a cy- 
clotron feature at —100 keV (|La Barbera et al.l2001[) giv- 
ing B ~ 10 13 G. Similar to SMC X-l, the soft component 
pulses more smoothly than the hard component and is 
perha ps out of phase with it l|Woo et alll996l |Pju i a_ej i jil] 
2002). For both SMC X-l and LMC X-4, iPaul et alJ 
(2002) argue that the soft component is likely not ther- 
mal in origin and is best described by power-law emission 
from the upper accretion column. 

XTE J011 1.2-7317 is a bright, transient XBP in the 
SMC with a probable Be star companion. It was dis- 
covered in 1998 Novemb er, simultan eously by RXTE 
(iChakrabartv et al.ll998|) and BATSE {Wilson fc Fingerl 
1998), and has a pulse period of 30.95 s. The source 
showed two outbursts of ~20 and 40 days in dura- 
tion, and the pulsar was spinning up with a very short 
timescale of P/P ~ 20 yr {Yokogawa et al J 12000). The 
likely Be companion m ay be embedded in an H II region 
{Coe et al.l2003 . 2000). Given a distance of 65 kpc to the 
SMC, the source had a 0.7-10 keV luminosity of -2 x 10 38 
ergs s _1 . The ASCA spectrum is an inverse, broken 
power law with iron emission at 6.4 keV, along with a 
soft component that was modeled as a b roken power law, 
therm al bremsstrahlung, or blackbody ( Yokogaw a et alJ 
2000). Variations in the soft X-ray band indicate that 
the soft component is pulsing. 

RX J0059.2-7138 is another bright transient Be/XBP 
in the SMC. It was discovere d as a serendip itous source in 
a short pointing of ROSAT {HughesHl994|) . The ROSAT 
plus ASCA spectrum has luminosity ~3 x 10 38 ergs s _1 , 
and the hard co mponent is fitted b y a cutoff power law 
with T - 0.4 {Kohno et al J 120001) . The soft compo- 
nent, for which L so f t /Lx ~ 35%, was fitted using ei- 
ther a soft broken power law in addition to the hard 
power law or a th in-thermal emission model (MEKAL; 
iMewe et al1ll986|) with fcT ME K = 0.37 keV. Fitting this 
model required a very low value for overall metal abun- 
dance in the emitting gas (Z = 0.015 Zq), but with a 
highly enhanced oxygen abundance in the absorbing col- 
umn {Zq ~ 7 Zq). Unlike the four sources discussed 
above, its soft component does not pulsate. 
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2.3. Other pulsars with a soft excess 

We know of three other well-studied XBPs that show 
a soft excess. 4U 1626-67, like Her X-l, is one of the 
few XBPs with a low-mass companion. X-ray pulsa- 
tions with period 7.7 s were first detected by SAS 3 
ijRappaport et alJll977|), and a faint, blue optic al coun- 
terpart was identified by |McC!mtoci^OD(Ez3)- There 
is an unusually short o rbital period of 41.4 minutes 
ijMiddleditch et aTlll981[K which rules out a hydrogen- 
burning main sequence star for the companion. In 1990 
June the steady spin-up of the source suddenly changed 
to a spin-down state that continues to the present. The 
spectrum has repeatedly been fitted with power law (r ~ 
0.7) plus blackbody (/cTrr ~ 0.3 keV) and emission-line 



features at ~1 keV 


Angclini et all 11995: 


Owens et alJ 


119971 lOrlandini et al. 


1998U ISchulz et alJ 


2001|). The 



pulse fraction is m uch lower for soft than hard X-rays 
ijOwens et al.111997)). suggesting that th e soft component 



is not pulsing. lOrlandini et al 
clotron absorption feature at 
B 



1998b) discovered a cy- 
35 keV, implying that 
3.5 x 10 12 G. They found the broadband (0.1-200 

34 



keV) X-ray luminosity to be (7.7 x 10"" ergs s 



- 1 ) D 2 



The distance is not well constrained, so it is difficult to 
make reliable conclusions about the physical processes in 
this source. 

There is also evidence for a soft excess in Cen X-3. It 
is a high-mass XPB, wit h pulse period of 4.8 s and or- 
bital period of 2.1 days ijSchreier et al.llT972b|) . and the 
O-ty pe companion is e stimated to be at a distance of ^8 
kpc (|Krzeminskilll974l) . The hard continuum has a typ- 
ical cutoff power-law shape, with iron line emissio n at 
~6.5 keV i|Nagase et alJll992 IBurderi et all 1200(3) and 
a cyc l otron line feature a t ~30 keV llSantangelo et al.l 
IT991 IBurderi et all 1200(1 ICoburn et alJ 120021) . giving 
B ~ 3 x 10 12 G. Using BeppoSAX data. IBurderi et al.l 
( 2000) found a bright pulsating soft excess with a black- 
body shape (A:Tbb ~ 0.1 keV). In this observation Cen 
X-3 had a large Lx ~ 10 38 ergs s _1 . 

Vela X-l also shows a soft excess. It h as a long pulse 
period of 283 s llMcClintock et al1 ll976f) and an orbital 
period of 8.96 davs llForman et alJH973f). The power-law 
spect rum has T ~ 1-1.7 l|Sako et all 119991 IChoi et al " 
1996) and the distance is ~1.9 kpc l)SadaklffleTt"ld 



indicating Lx ~ 10 ergs s . There are cy- 
clotron features at ~ 24 and 5 5 keV l)Krevkenbohm et al.l 
l2002HKretschmar et aLlll996j) . indicating B ~ 2.4 x 10 12 
G. The ROSAT spectrum showed a soft, non-pulsing 
thermal bremsstrahlung component with kT t ^ ~ 0.5 keV, 
attributed to shock -heated optic ally thin emission by the 
strong stellar wind (Haberl 1994). This was similar to the 
high-mass X-ray binary 4U 1700-37, which had a remark- 
ably similar R OS A T soft excess spectrum l)Haberl et al.l 
1994). However, high-resolution grating spectroscopy 
of these sources with Chandra showed that in eclipse 
there were emission lines that likely came from material 
that was not collisionally excited bu t was photoionized 
by the X-rays from the neutron star ijSchulz et al.ll2002T: 
IBoroson et al.ll2003|) . We discuss collisionally heated and 
photoionized gas emission in §§ 5-6, respectively. 

There are two additional transient Be/XBPs in the 
LMC for which a spectrally-fitted soft excess has been 
seen. EXO 053109-660 9.2 was the b r ightes t source in 
a survey of the LMC bv lHaberl et all (J2003), who esti- 



mated Lx — 4.6 x 10 ergs s _1 . They fitted the soft com- 
ponent as thin-thermal (MEKAL) emission with reduced 
abundances (0.29 Z@) of elements heavier than oxy- 
gen. A 0538-66 was seen in two outbursts with ROSAT 
ijMavromatakis fc Haberll Il993). when the source had 
Lx ~ 2-4 x 10 37 ergs s _1 . They fitted the spectra with 
a power law, T ~ 0.8, a soft blackbody, kT^ ~ 0.2 keV, 
and a thermal bremsstrahlung, kT t b ~ 0.25, although 
the spectra have a relatively low number of counts so 
there may be considerable uncertainty in these parame- 
ters. The authors speculate that the blackbody compo- 
nent may be emitted from near the Alfven radius of the 
neutron star, while the bremsstrahlung component may 
be powered by collisional heating of the stellar wind. Al- 
though we mention EXO 053109-6609.2 and A 0538-66 
here, they have little spectral information in soft X-rays, 
so we do not discuss them in detail. 



3. UBIQUITY OF THE SOFT EXCESS 

In this section we address whether the soft excess is 
a universal feature of XBPs. All the brightest sources 
with low absorption show this feature, so it is natural to 
suspect that it may be present also in other XBPs but 
not be detected because of either low flux or high absorb- 
ing column. In Table 2 we show spectrally-determined 
Ah values and the unabsorbed fluxes for a number of 
known XBPs. The sample is taken from Table 1 in 
IBildsten et all ltl997l) . plus EXO 053109-6609.2. We 
have included only sources for which there are observa- 
tions with sensitivity at the energies of the soft excess 
less than 2 keV. The fluxes are defined in different bands 
in the range 0.1-20 keV, so there is some error (factors of 
a few or so) in the relative fluxes from these sources, al- 
though this does not affect the overall pattern observed. 

The sources in the sample are plotted as a function of 
flux and Ah in Fig. 1. The XBPs with a known soft 
excess are shown as open stars; those without as solid 
squares. It is clear that the brighter systems with less 
absorption are far more likely to show a soft excess, and 
that no sources with high absorption (Ah > 3 x 10 22 
cm" 2 ) show this feature. There is an apparent boundary 
between the sources that show a soft excess and those 
that do not, suggesting that the observability of this fea- 
ture depends on flux and Ah- 

Marked with an asterisk is the Be/XBP X Persei, 
which ha s a markedly d i fferent soft excess from the other 
sources. ICoburn et alJ l)2001|) modeled this component 
as a blackbody with &Tbb ~ 1.4 keV, which is an or- 
der of magnitude higher than the typical temperature. 
X Per is very nearby ( D ~ 0.7 kpc) and un der- luminous 
(L x ~ 10 34 ergs s" 1 ) l|di Salvo et alJH998j) . and the au- 
thors attributed the soft component to blackbody emis- 
sion from the accreting polar cap. 

We modeled the detectability of the soft excess at 
various fluxes and Ah by performing fits to simulated 
data. The response of the ROSAT PSPC detector 
(Pfcff ermann fc Bri el 1986) is used, because many XBPs 
have been observed with this detector and beca use it is 
sensit ive in the soft X-ray band. Using XSPEC ijArnaudl 
1996), we created simulated PSPC spectra of a source 
with a spectrum like that of SMC X-l (r = 0.9, kT B B = 
0.15 keV, and L so f t /Lx = 0.04), and varied Ah and the 
overall unabsorbed flux in the 0.2-3 keV band. The sim- 
ulated exposure time was 25 ks, and counting statistics 
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TABLE 2 

Unabsorbed flux AND TVh values for a sample of XBPs 



Source 


Flux (ergs cm 2 s 1 ) Energy range (keV) a 


N H (10 22 cm" 2 


) Mission Reference 




Soft excess observed 








Her X-l 


3.7 X 10 -8 0.3-10 


0.005 


XMM 


1 


SMC X-l 


4.7 x 10~ 10 0.7-10 


0.5 


ASCA 


2 


LMC X-4 


4.0 x 10 -10 0.7-10 


0.057 


ASCA 




XTE 0111.2-7317 


3.6 x 10~ 10 0.7-10 


0.18 


ASCA 


4 


RX J0059. 2-7138 


6.9 x 10 -11 0.2-2 


0.088 


ROSAT 


5 


4U 1626-67 


6.4 x 10~ 10 0.5-10 


0.5 


Chandra 


6 


Cen X-3 


5.7 x 10 -9 2.0-10 


1.95 


BeppoSAX 


7 


Vela X-l 


5.1 x 10~ 9 2.0-10 


0.86 


BeppoSAX 


8 


EXO 053109-6609.2 


1.5 x 10 -10 0.2-10 


0.69 


XMM 

J\ Ivl Ivl 


9 


A 0538-66 


1.3 x 10~ 10 0.1-2.4 


0.08 


ROSAT 


10 




So ft GXCGSS Tiot obSGTVGO 








riRO 11 744— 98 


9 3 V 1 — 8 9 D— 1 n 


5 


ASCA 


11 


GX 1+4 


8.4 x 10" 10 2.0-20 


20 


ASCA 


12 


OAO 1657-415 


4.7 x 10~ 10 2.0-10 
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Fig. 1. — Observed A^h an d unabsorbed fluxes for the XBPs in 
Table 1 of lBildsten eTantT99^ . plus EXO 053109-6609.2. Sources 
with a known soft excess are shown as stars, while those without 
are shown as squares. The low-luminosity source X Per is shown 
as an asterisk (see text for details). 



Fig. 2. — Observability of the soft excess as a function of Nn and 
unabsorbed fluxes, from fits to simulated ROSAT data. Contours 
of x% f° r n t s to a simple absorbed power law are shown. The soft 
excess is observable for high \i an d n °t observable for low xt- The 
XBPs from Fig. 1 are also shown for comparison; note that the x 2 
contours are similar in shape to the boundary between observed 
XBPs with and without a soft excess. 



were included. We also included a typical background 
spectrum, taken from the PSPC observation of SMC X- 
1 from 1991 October. This background has a surface 
brightness in the PSPC band (0.1-2.4 keV) of 2.1 x 1(T 2 
counts s^ 1 arcmin -2 , which is similar to those found for 
Cen X-3, LMC X-4, and Her X-l. More details on PSPC 
observations of these sources are given in § 4.1. 
Next, we determined whether each spectrum, after 



background subtraction, could be fitted with an absorbed 
power law without a soft excess component. For the fits 
we took initial values of Nn and T and overall normal- 
ization from the simulated data, and allowed all these 
parameters to vary. We recorded \t f° r the best fit to 
each simulation. For low xt> we conclude that a power 
law spectrum satisfactorily fits the simulated data and 
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Fig. 3. — Schematic showing the four soft excess emission processes considered in this paper. 



so the soft excess is not detectable. For high xt the soft 
excess is detectable. 

The fit results are plotted in Fig. 2, along with the 
XBPs from Fig. 1. Contours of xl = 1-2, 1-5, and 2 
are shown. For low flux and high Ah (lower right), xt 
varies somewhat randomly between 0.8 and 1.2; for clar- 
ity, these contours are not shown. We find that the soft 
excess is detectable in sources with high flux and low 
Ah, and undetectable for low flux and high Ah- The 
steep contours of Xui which indicate the boundary of de- 
tectability, are close to and similar in shape to the actual 
boundary seen in observations. This suggests that many 
sources with low flux or high Ah do in fact have a soft 
component but that this feature has not been observed 
because of selection effects. We conclude that the soft 
excess is probably a very common intrinsic feature of 
XBPs. 

As discussed in § 2, observations of the soft excess in 
XBPs have been modeled in a variety of ways, many 
of which require different emission processes. Given the 
similarities of many of the observed soft components and 
the apparent ubiquity of this feature, it is likely that it 
has a common origin in many sources. Possible mecha- 
nisms can be separated into four main types: 

1. emission from the accretion column, 

2. thermal emission from collisionally energized, dif- 
fuse gas around the neutron star, 

3. reprocessing of the hard (power law) X-rays by dif- 
fuse material around the neutron star, and 

4. reprocessing by optically thick, dense material 
around the neutron star. 

In the following sections, we consider these four scenar- 
ios and examine them in light of observational data and 
the physics of the relevant emission. Simple schematics 
of the various models are shown in Fig. 3. For a model 
to be viable it must satisfy two constraints: 



1. the energy balance of the system must be able to 
produce the observed soft luminosity and 

2. the spectrum, pulsation properties, and other ob- 
served features must be consistent with model pre- 
dictions. 

4. EMISSION FROM THE ACCRETION COLUMN 

It is natural to suspect that the soft component origi- 
nates in the same region as the hard component, that is, 
in the column of accreting gas. Low-lumin osity soft emis- 
sion f rom the polar cap may exist in X Per l)di Salvo et alJ 
1998), and the more luminous soft components in SMC 
X-l, LMC X-4, and XTE J0111. 2-7317 have been mod- 
eled by non-thermal broken power-law em i ssion from 
the a ccretion column ( Yokoeaw a et al.ll2f)0ft iPaul et all 
|2002|) . This model predicts soft pulses, because emis- 
sion is confined to a small region that rotates with the 
neutron star. 

4.1. Brightness temperatures of accretion columns 

In XBPs the accreting material flows along magnetic 
field lines from the magnetosphere, at radius R m , to the 
surface of the neutron star at -Rns- To estimate the 
surface area of the accretion column, we follow § 6.3 of 
iFrank et al.l (|2p02). We assume that the magnetic field 
is a dipole that is tilted by an angle a with respect to 
the accretion disk, and that only field lines that pass 
through the disk outside R m accrete. It is straightfor- 
ward to show that all the accreting field lines intersect 
the neutron star inside an angle (3 from the pole, where 

sin 2 f3 = (RNS/Rm) sin 2 a. 
Thus the radius of the accreting polar cap is 

i? co i — i?NS sin/3 ~ -Rns(-Rns / Rm) 1 ^ 2 , 

assuming sin a ~ 1. For the XBPs we consider, R m ~ 
10 8 cm so for i?NS ~ 10 6 cm, i? co i ~ 0.1i?Ns- If we 
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TABLE 3 

ROSAT PSPC DATA USED FOR Tf, CALCULATIONS 



Source 


Date/Duration (ks) 


(10 22 cm" 2 ) 


Her X-l 


1993 Aug 


5 


0.003 


SMC X-l 


1991 Oct 


17 


0.36 


LMC X-4 


1991 Oct-Nov 


33 


0.05 


RX J0059. 2-7138 


1993 May 


5 


0.04 


Cen X-3 


1992 Jan 


9 


2 


Vela X-l 


1992 Nov 


17 


0.5 



picture the accretion column as a cylinder with height 
< i?NS ) the surface area of the column is 

Acoi 27ri? co ii?NS = /coi47ri?^, s . 

Here / co i is simply the surface area of the accretion col- 
umn as a fraction of the neutron star's surface area, with 
/col < 0.1. 

For high L S oft, this small area implies a very high tem- 
perature for the emitting gas. The minimum tempera- 
ture for material radiating a spectral luminosity L„ in 
a given band is the brightness temperature Tf,, which is 
the temperature of a blackbody that would produce the 
same L v . For an emitting area A, 

L v = AnB u (T b ). 

The value of Tf, for different energy ranges sets a lower 
limit on the temperature of the gas emitting at those 
energies. 

We used data from the ROSAT PSPC observations for 
Her X-l, SMC X-l LMC X-4, and Cen X-3, as well as RX 
J0059. 2-7138 and Vela X-l, which we include for com- 
parison even though their lack of soft pulsations would 
appear to rule out accretion column emission. Data on 
these observations are given in Table 3. The PSPC is a 
proportional counter detector with response at low ener- 
gies, from 0.1 to 2.4 keV, and wit h spectral resolution 
AE/E = 0.43(£:/0.93 keV)" 1 / 2 fPfefferma rm fc Briell 



1986). We obtained the data from the NASA HEASARC 
archive and used the standard event-screening criteria. 
We extracted the source and background spectra from 
circular and annular regions centered on the source. In 
the case of RX J0059. 2-7138, a nearby supersoft source 
(IE 0056.8-7154) wa s excluded from th e selection re- 
gions, as described in lKohno et alJ 1)20001 ). 

We assume that all the observed flux in the spectrum 
comes from a column of area j4 co i and determine Tf, as a 
function of energy. We use energy bins corresponding to 
the PSPC resolution and fit each bin with a power law 
of constant energy (photon T = 1), modified by typical 
values of V H (shown in Table 3). Each XSPEC fit gives 
the normalization of the power-law I p \ (in photons cm~ 2 
s _1 keV -1 at 1 keV) and its error. The energy flux in 
each bin (from E\ to E^) is then 

F ohs = / E(I pl E- 1 )dE = I pl (E 2 -E 1 ), 

J Ei 

Given distances listed in Table 1, we calculate the lumi- 
nosities and then Tf, in each bin. To obtain a lower limit 
on T b , we use the 3<r lower limit on 2pl. 
The isotropic luminosity observed in a given bin is 

Lobs = 47rD 2 F obs , 
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" RX J00h9 _ 
" Cen X-3 







0.0 0.5 1.0 1.5 2.0 
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Fig. 4. — Brightness temperatures for the observed fluxes (3ct 
lower limit) from ROSAT PSPC data for a sample of XBPs. We 
assume an emission region with area a fraction / GO ] = 0.1 of the 
surface area of the neutron star (Rns = 10 km). For E <g Tf,, 



and the luminosity emitted from the accretion column in 
that bin is 



A col ir / B(T b )dE. 

J E x 

Letting L obs = L cm , and using A co i = / co i47rii|, s , we 
have 

D 2 

B(T b )dE — — — g — - — -Fobs- 

We solve this numerically to find Tf,. In the Rayleigh- 
Jeans limit E <C kT b this gives, for / co ] = 0.1, i?NS = 10 
km, Ei and E% in keV, and L38, the luminosity in the 
bin in 10 38 ergs s _1 , 



T b = 1.75 x 10 10 K 



L 



El 



Ef 



For each source we plot the values of T b , given f co \ = 0.1, 
in Fig. 4. Vela X-l and Cen X-3 are highly absorbed be- 
low 0.5 keV, so we do not calculate Tf, for those energies. 

As expected, Tf, is higher for more luminous sources. 
For low energies, we find that for SMC X-l, LMC X-4, 
and RX J0059. 2-7138, the observed flux gives T b > 10 11 
K, even for a relatively large accretion column area 
(/coi = 0.1). At higher energies, we have Tf, > 10 8 K 
(and, remarkably, T b > 10 9 K for SMC X-l). These T b 
values represent the minimum possible temperatures of 
the gas emitting at these energies. Are these temper- 
atures consistent with models of accretion columns? In 
general, such models have been separated into two types, 
depending on the shock mechanism that decelerates the 
gas as it approaches the neutron star's surface. The 
shock can be dominated by gas pressure (for low Lx sys- 
tems) or by radiation pressure (for Lx > 10 36 ergs s _1 ). 
The latter models are applicable for the luminous XBPs. 
This type of accretion column consists of a mound of 
subsonic, set tling materia l below a radiation-dom inated 
shock front l)Beckerlll998(l . IBurnard et"all 1)1991(1 mod- 
eled the dynamics and emission properties of the accre- 
tion mounds, and found that temperatures in the mound 
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Fig. 5. — Spectra of c yclo-synchrotron r a diatio n for different 
temperatures, taken from [Mahadcvan ct afj 119961) . Spectra are 
plotted as a function of x = E/Eg. 



vary from 10 s K at the outer boundary to 2 x 10 8 K at 
l arge optical dep t hs. 

Bur nard et alJ l)1991fl note that the accretion mound 
emission is Compton-scattered by the overlying shock, 
which could increase 71,. In principle, if the gas free-falls 
all the way to the neutron star's surface and is stopped 
in a t hin shock, then temperatur es of ~10 12 K are pos- 
sible dShapiro fc Teukolskvlll983(l . However, we expect 
that collisions in the column would slow the gas, and it 
is expected that a radiative shock is not perfectly thin 
but has w idth approxim ately equal to the radius of the 
polar cap ( Becker 1998). Therefore we expect the shock 
temperature to be lower. It is difficult to see how Comp- 
tonization by a real shock could produce the very high 
soft fluxes observed. 

This result raises the interesting, puzzling question of 
how, for SMC X-l, even the hard emission above 1.5 keV 
could be produced by the accretion column. For a large 
column size / co i = 0.1, the Tf, ~ lfj 9 K is hotter than 
models predict. This suggests that either the column is 
larger in extent (i.e. f co \ > 0.1) or it is more energetic 
(T > 10 9 K) than predicted. This question is worth 
investigating in future work. 

4.2. Comparison with cyclo-synchrotron spectra 

Even if accretion column temperatures are very high, 
it is difficult to produce the observed soft spectral 
shapes. For high Tf, > 10 7 K, blackbody or thermal 
bremsstrahlung emission would produce the majority of 
the flux at energies ~ kT > 1 keV, which is incon- 
sistent with observations. One possible way out is to 
note that, in the large magnetic field of the accretion 
column, the dominant emission process is likely to be 
cyclo-synchrotron radiation. If this radiation is thermal, 
the peak of the emission can be at much lower energies 
than kT. The shape of this spectrum has been calculated 
in rigorous detail bv lMahadevan et alJ |)199(tT ). who show 
that it is a function of x — E/Eb, where 

Eb = ~J~~~ = H.6(S/10 12 G) keV 

2nm e c 



is the (non-redshifted) cyclotron energy. Therefore, 
the energy range of cyclo-synchrotron emission depends 
largely on B. We plo t the calculated spectra from 
iMahadevan et all l)1996j) in Fig. 5. 

AtT = 3.2xl0 10 K, which is smaller than the mini- 
mum Tb that we obtained for the luminous XBPs at low 
energies, the spectrum peaks at logy ~ 3. A typical soft 
component in the luminous XBPs has the peak of vF v 
at ~ AkT-QB ~ 0.7 keV. For cyclo-synchrotron emission 
to peak at this energy, we require that 

E B ~ 0.7 keV/x - 7 x 10~ 4 keV, 

implying that B ~ 6 x 10 7 G. For a dipole field with 
surface strength 4 x 10 12 G, 

B(r) ~ 4 x 10 12 {r/R m r 3 , 

implying that the emission must come from a radius r ~ 
40i?NS- This is much larger than the expected size of the 
accretion column, and the gas cannot attain such a high 
temperature so far from the neutron star. 

It is difficult to conceive of any other emission mech- 
anism that could produce a high soft luminosity from a 
small area and yet have the spectrum cut off at the ener- 
gies observed. We conclude that in the luminous XBPs 
(SMC X-l, LMC X-4, RX J0059.2-7138, XTE J0111.2- 
7317, and Cen X-3), the soft component must be emitted 
from a region larger than the surface of the neutron star. 

5. EMISSION BY A COLLISIONALLY ENERGIZED CLOUD 

A large, diffuse cloud of gas around the neutron star 
would have the size needed to produce the high L so ft 
and might be powered by collisional energ y. Such emis- 
sion may exist in Vela X-l a nd 4U 1700-37 (Schul z~t alJ 
120021 iBoroson et alJ 12003), and perhaps also in RX 
J0101.3-7211 and AX J0103-722, which were observe d 
in the XMM survey of the SMC bv lSasaki et all (|2003|) . 
These sources have Lx ~ 10 35 ergs s -1 , and were fitted 
with a power law plus an optically thin thermal emission 
model (MEKAL) with kT MEK ~ 0.2 keV. 

Collisional heating is expected in XBP systems; a 
detailed hydr o dyna mic simulation of LMC X-4 by 
IBoroson et alJ l)2001[) found shocks and wake structures 
that coul d heat the gas and cause optically thin X-ray 
emission l)Fransson fc Fabianlll980f ). Fitting a thin ther- 
mal emission model to data from SMC X-l and LMC 
X-4, iPaul et all {2002) found that the observed high 
L S oit requires n 2 V ~ 10 61 , which for typical gas den- 
sity n < 10 12 cm -3 implies a very large emission region 
of -Rcioud > 10 12 cm. The same is true for XTE J0111.2- 
7317 and RX J0059. 2-7138, which have similar soft spec- 
tra. 

The large emission region immediately rules out this 
process for sources with a soft component with short 
pulsations, such as SMC X-l and XTE J0111. 2-7317. 
A pulse period of about a few seconds is much shorter 
than the light-crossing time of such a large cloud, and 
so the cloud's emission could not pulsate. For SMC X- 
1 and LMC X-4, the size of the emitting region is also 
constrained by the drop of the soft flux during eclipses 
ijVrtilek et alJ I200H IBoroson et al1l2001|) . The eclipse 
light curves for SMC X-l and LMC X-4 fe.g.. lWoo et all 
ll995HLa Barbera et al.l200lD show that the soft flux falls 
with a timescale < 10% of the total eclipse duration. 
Thus, i?cioud < 0.1i?+ < 10 11 cm. Such arguments rule 
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out collisionally energized emission for SMC X-l, LMC 
X-4, and XTE J011 1.2-7317, but the process may still 
be important in RX J0059. 2-7138, which shows neither 
eclipses nor pulsation at soft energies. 

In XBPs the main input of kinetic energy is from a 
wind, simply E — Mv^/2, where v w is the wind veloc- 
ity, typically ~2000 km s _1 . Mass loss rates from the 
star and from the irradiated accretion disk are typically 
M < 1CT 6 M Q yr-\ giving E < 5 x 10 36 ergs s -1 . Such 
emission can easily power the soft components of lower 
luminosity XBPs such as Vela X-l. However, it is possi- 
ble for the more luminous soft excess of RX J0059. 2-7138 
only if there is a very vigorous wind for which all the ki- 
netic energy is converted to soft X-rays. This is unlikely. 

In addition, diffuse thermal emission in RX J0059.2- 
7138 should show line emission fro m metals. In their 
MEKAL fit to the ASCA spectrum. Ikohno et all (2000) 
obtain an extremely low Z — 0.02, which fits best be- 
cause it greatly reduces the flux of most X-ray lines. 
While the metallicity of the SMC is expected to be low, 
this value is unrealistic. It is more likely that the non- 
detection of lines is due to optically thick (blackbody- 
like) emission, which we discuss in § 7. 

We conclude that collisionally energized diffuse emis- 
sion can power the soft excess only in systems with 
^ 10 36 ergs s _1 . This may explain part of the 
soft components in Vela X-l and 4U 1700-37, but in 
RX J0059. 2-7138 the soft excess is too luminous for this 
process. With high spectral resolution observations, it 
is possible to easily identify this type of emission by the 
presence of X-ray emission lines. In addition there should 
be other spectral signatures from diffuse gas around 
XBPs, because of photoionization of the gas by the cen- 
tral X-ray source. We discuss this process in the next 
section. 

6. REPROCESSING BY A DIFFUSE CLOUD 

It is possible that the soft flux arises from reprocessing 
of the hard X-rays by a diffuse cloud of gas around the 
neutron star. Such a cloud must absorb and re-radiate a 
fraction of the incident hard luminosity equal to / rep = 
Lsoft/Lx ~ 0.1 and must produce an output spectrum 
like those observed in XBPs. Because of the high Lx in 
these systems, the gas is highly photoionized by the X- 
ray flux. Emission from such a gas may co ntribute to the 
soft excess in Vela X-l and 4U 1700-37 (iBoroson et al.l 
l2003t iSchulz et al.ll2002]) . The XSTAR code is designed 
to calculate the detailed structure of photoionized gas, 
and has been used extensively to model the circumstellar 
environments in XBPs. For details on the c ode and its 
applications, sec Kallm an fc Bautistal l)2001|h 

XSTAR works by determining the temperature, ioniza- 
tion and excitation structure of a spherical gas cloud illu- 
minated by a central point source of X-rays. It calculates 
the absorption by the gas due to photoionization and 
Compton heating and the emission from the gas through 
recombination, free-free emission, and atomic line emis- 
sion. In the case of an optically thin low-density (< 10 12 
cm ) gas, the structure is a function simply of the ion- 
ization parameter £ = L/nR 2 , where L is the luminosity 
of the central source, n is the numbe r density of the gas , 
and R is the distance from t he source dTarter e t al.ll 969*1. 

XSTAR has been u sed bv lWoidowski et all ( 2000) and 
IBoroson et all i|2001f) to compare the observed eclipse 




Fig. 6. — Schematic of the XSTAR models. The observed emis- 
sion is taken to be the sum of the input power law and the outward 
emitted spectrum. 



spectra from SMC X-l and LMC X -4 with hydrodynamic 
simulations. iVrtilek et"aT1 l)2001|) used XSTAR to ana- 
lyze relatively low-luminosity (^10 35 ergs s _1 ) photoion- 
ized line emission from SMC X-l seen in Chandra ACIS 
spectra. In eclipse and in the low state (when the neu- 
tron star was blocked by the disk), they found lines that 
correspond to XSTAR simulations with 1 < log£ < 1.5. 
IVrtilek et aT] (|2001^ found that at a radius of 2 x 10 12 
cm from the neutron star, these values for £ indicate a 
gas density of ~10 12 cm -3 . 

For most of the XBPs with Lx > 10 38 ergs s _1 , we 
do not expect a diffuse cloud to be the source of the soft 
excess as discussed in § 5. However, such a scenario is 
possible in RX J0059. 2-7138, which does not show soft 
pulsations or eclipses, and so we test the possiblity that 
optically thin reprocessing is the main source of the soft 
excess in this system. We have used the latest version of 
XSTAR, version 2.1k, to simulate the photoionized gas 
in this source. 

6.1. XSTAR input 

Because our models must reprocess a substantial frac- 
tion of the incident X-ray flux, we required that the inter- 
vening gas subtend a large solid angle around the central 
source. We thus modeled the reprocessing by spherical 
shells of gas centered on the neutron star. For a cloud of 
radial extent less than 10 11 cm, column densities higher 
than 10 23 cm -2 imply densities higher than 10 12 cm -3 , 
for which we cannot parameterize our models simply in 
terms of £. 

Instead, we constructed a series of models of gas in 
spherical shells with varying radius R and column density 
N (Fig. 6). Because the gas density in the system falls 
off with radius, the largest clouds of constant density 
should have thickness comparable to the scale length, or 
AR ~ R. We therefore ran models with AR — R and 
OAR to explore the effects of higher density parcels. 

For the input spectrum we used an exponentially cut- 
off power law that matc hes the hard c ompo nent of RX 
J0059. 2-7138 as fitted bv lKohno et"afl pOOOh . The spec- 
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TABLE 4 
Results of selected XSTAR RUNS a 



logB 


logN 


log n 


bj 






z 

;T 


= 0.25 

/rc 

rep 


z 

logT 


= 0.01 

/rc 

rep_ 




z = 

l°gT 


1 

/rc 
rep 
















AR = R 












9.0 


24.0 


15.0 


4.7- 


-5.4 


7.6- 


-7.6 


1.0 x 10" 4 


7.6-7.6 


1.3 X 10" 


-5 


7.6-7.6 3. 


8 X 10" 4 


10.0 


24.0 


14.0 


3.7- 


-4.4 


6.9- 


-7.5 


3.6 x 10" 3 


7.0-7.5 


4.9 X 10" 


-4 


6.9-7.5 


0.014 


11.0 


23.5 


12.5 


3.2- 


-3.9 


4.6- 


-7.2 


0.045 


4.6-7.2 


0.099 




6.4-7.2 


0.018 


11.0 


24.0 


13.0 


2.7- 


-3.4 


4.7- 


-6.5 


0.14 


4.5-5.6 


0.022 




6.2-6.6 


0.14 


11.5 


23.5 


12.0 


2.7- 


-3.4 


4.6- 


-6.4 


0.087 


4.3-5.1 


0.040 




5.2-6.6 


0.099 


12.0 


23.0 


11.0 


2.7- 


-3.4 


5.4- 


-6.5 


8.7 x 10" 3 


4.8 b 


7.7 x 10" 


-3 


6.2-6.6 


0.022 
















AR = 0.1R 












11.0 


24.0 


14.0 


2.3- 


-2.4 


4.5- 


-5.1 


0.21 


4.5-4.9 


0.045 




4.4-5.8 


0.41 


11.5 


23.5 


13.0 


2.3- 


-2.4 


4.1 


i h 


0.14 


4.7 b 


0.044 




4.6 b 


0.26 


12.0 


23.0 


13.0 


2.3- 


-2.4 


4.1 


S b 


0.061 


4.5 b 


0.051 




4.6 b 


0.14 



a All quantities are in cgs units. 

b Modcl run with constant temperature throughout cloud. 



trum has luminosity Lx — 2.6 x 10 38 ergs s -1 , T = 0.43, 
£cut = 6.4 keV, and -Efoid = 9.3 keV. For gas co mposi- 
tion we used a standard SMC value of Z = 0.25 ijWeltvl 
1200 ID . as well as Z = 0.01 and 1. 

There are large turbulent velocities in the gas around 
the neutron star, so we modified the standard XSTAR 
code to include these effects in the line opacity. For each 
model we set the turbulent velocities equal to the Kep- 
lerian velocity, 

wturb = v K = 1.9 x 10 4 (i?/10 8 cmT 1/2 km s _1 
for a neutron star mass of 1.4 M Q . This decreases the 
line opacities by a factor of fthermai/(^turbV A A) where A 
is the atomic mass. 

We ran models with radii from 10 s cm (the expected 
radius of the inner disk) to 10 12 cm (a few times the size 
of the binary system during outburst which is R ~ i?* ~ 
5Rq). Most of the reprocessing column should be con- 
tained within this radius. Close to the neutron star, the 
accretion flow can create large gas densities and columns 
^10 24 cm -2 or greater. At R > the wind from the 
Be star is expected to blow away the surrounding in- 
terstellar medium so that the gas density is determined 
mainly by the star's mass loss. For a spherical wind, 

M = n/XTOH47ri? 2 w w 

where fj, is the molecular weight, ~1.4. For M < 10~ 6 
M Q yr _1 and v w > 1000 km s _1 , we have upper limits in 
our models of N < 10 23 5 and 10 23 cm" 2 , for R = 10 115 
and 10 12 cm, respectively. We note that while large den- 
sities may be possible at smaller radii, we did not run 
models with N > 10 24 cm" 2 because the effects of Comp- 
ton scattering are not explicitly included in XSTAR. We 
discuss Compton-thick reprocessing in the next section. 

In most models, we allowed XSTAR to vary the tem- 
perature in each subshell to minimize heating minus cool- 
ing as it stepped outward through the cloud. However, 
for some runs at large radii (10 115 — 10 12 cm), the model 
encountered problems converging on the temperature. In 
these cases we set the entire cloud to be at the single tem- 
perature that minimized heating minus cooling over the 
entire cloud. 



6.2. XSTAR results 

Outputs of sample model runs are given in Table 4. 
Because a simple power law fits the hard end of the spec- 
trum well, there is likely to be little absorbing gas in the 
system directly along the line of sight, so the observed 
spectrum is the sum of the input hard power law plus 
the reprocessed emission from surrounding gas. Samples 
of these spectra are shown in Fig. 7. 

We find that none of the models are able to reproduce 
the observed emission from RX J0059. 2-7138. In models 
in which log£ > 5 (that is, for small radii) there is vir- 
tually no reprocessing for any column density, because 
the gas temperature is very high (~5 x 10 7 K), all the 
atoms in the gas are completely ionized, and therefore, 
the opacity to photoionziation interactions is low. Repro- 
cessing inside this radius requires densities much higher 
than those that are possible to model using XSTAR. This 
approaches dense, optically thick reprocessing, which is 
discussed in the next section. 

For larger R > 10 11 cm (and thus smaller £), the 
spectrum is markedly altered by the photoionized gas, 
as shown in Fig. 7. Some of the hard X-rays are 
absorbed through ionization and re-emitted as softer 
photons through recombination, bremsstrahlung, or line 
emission. For each model, we calculate / re p) the ratio 
of the total luminosity emitted outward by the gas to 
the total input luminosity. For most models / rep < 0.1, 
and at most 50% of the reprocessed emission emerges in 
soft X-rays (0.1-2 keV). There are a few cases for the 
AR = OAR shells that have / rop > 0.2, but in these 
only 10%-20% of the emission emerges in soft X-rays, 
with most coming in the EUV. In addition, the spectral 
shape of the soft emission does not match observations. 
There is little bremsstrahlung continuum, and the ra- 
diative recombination and atomic lines produce spectra 
with discrete features and insufficient flux at ~1 keV. 

To show this quantitatively, we used the output XS- 
TAR spectra to create table models in XSPEC and fitted 
these models to the observed spectrum of RX J0059.2- 
7138. We too k data from the A SCA GIS2 observation 
of 1993 May l|Kohno et all l20"o7l from the HEASARC 
archive. We used the standard event-screening criteria 
and extracted the source and background spectra from 
regions similar to those shown in Fig. 1 of iKohno et alJ 



SOFT EXCESS IN X-RAY PULSARS 



11 




Fig. 7. — Left: Sample photon spectra from XSTAR showing the effects of increased reprocessing with larger radius. The models have 
Z = 0.25 Zq, N = 10 cm -2 , and AR, = R. The dotted line shows the input spectrum, matching the exponentially cutoff power law for 
RX J0059. 2-7138. Note the recombination continuum at energies below 0.1 keV. Right: Part of each spectrum for 0.7-10 keV, convolved 
with the response of the ASCA GIS2 detector. Note that features become visible in the ASCA spectrum as reprocessing increases. 



( 2000). We fitted the data in the range 0.7-7 keV. 

For each XSTAR run, the model consisted of the input 
power law spectrum plus the outward emitted spectrum, 
modified by typical interstellar Ah = 4 x 10 20 cm~ 2 . 
We fitted these models to the data, allowing the rela- 
tive luminosities of the input and emitted component to 
vary between 0.01 and 100 times that calculated by XS- 
TAR (indicating more or less reprocessing gas than in 
our model shells). The models fitted, all having xt > 2.5 
and most having xt > 7- In general, the reprocessed 
spectra have discrete line features and too little flux at 
~1 keV. In addition, most spectra are underluminous; 
for the model with the best fit (Z = 0.01, R = 10 12 
cm, N = 10 23 cm -2 , and AR = 0.1R), the condition 
xt = 2.5 requires an outward luminosity 90 times that 
output by XSTAR. 

To allow for reprocessing by gas at several different 
radii and densities, we also performed spectral fits with 
multiple models simultaneously. For each Z we included 
the emitted spectra from the six models with the most 
luminous emission in soft X-rays (0.1-2 keV) plus the 
input power law. We allowed the normalization of each 
emitted spectrum and the input power law to vary, and 
the results are shown in Fig. 8. Even when combin- 
ing emission from several sets of gas parameters, there 
remains a decrement at low energies, which is not signif- 
icantly affected by changing Ah- As shown in § 4, this 
required soft flux could not simply be part of the emission 



from the neutron star because of the high temperatures 
required. 

Of course, the photoionized gas in XBP systems is 
far more complex than we have modeled here, contain- 
ing shocks, wakes, and other density gradients, and in 
the case of Be stars, a circumstellar disk. We have not 
considered small regions of very dense gas (n > 10 12 
cm ) at large ra dii, such as those pre sent in simula- 
tions of SMC X-l llBlondin fc Wodll99ll) and suggested 
by the lines seen bv lVrtilek et all l|2001jK We note, how- 
ever, that such dense regions subtend a small angle to 
the neutron star and so do not contribute significantly 
to the reprocessed flux. There also could be variations 
in the shape of the photoionizing spectrum seen by dif- 
ferent parcels of gas, which would affect the photoion- 
ization emission in complex ways that would require a 
much more detailed study to model completely. Thus, 
it is in principle possible that a particular configuration 
might produce the soft spectrum seen in RX J0059.2- 
7138. However, our analysis suggests that in general 
photoionized emission is too weak and would have an 
incorrect spectral shape. In RX J0059. 2-7138, as in the 
other luminous XBPs, the soft excess most likely has a 
different origin. 

6.3. Higher optical depths 

Because XSTAR does not include Comptonization, 
our models have not covered the case of Compton-thick 
(N > 10 24 cm~ 2 ) gas at radii close to the neutron star 
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Fig. 8.— Observed ASCA GIS2 spectrum for RX J0059.2-7138, 
fitted with the input power law plus the six most luminous outward 
emission spectra from XSTAR, for Z = 0.25, 0.01, and 1 Zq. Each 
fit had 24 degrees of freedom, and the best-fit xt values are as 
shown. The dashed lines indicate the contributions from the power 
law and the reprocessed emission. Note that even when emission 
from multiple models is included, the spectra do not reproduce the 
observations. 



where high densities may be present. IRossI ()1979j) models 
this regime with code that explicitly calculates the ra- 
diative transfer in Compton-thick gas. The models show 
that for R ~ lfj 8 cm, high densities (n ~ 10 18 cm' 3 ), 
and Lx = 2 x 10 37 ergs s -1 , there is still some decre- 
ment from a blackbody shape at X-ray energies, even for 
high Thomson optical depths tt ~ 3-5. In this regime 
there are strong emission features of oxygen and silicon 
at ~ 0.8- 1.2 keV, which are not seen in the soft excess 
of most XBPs but may correspond to the ~1 keV emis- 
sion lines found in Her X-l. This suggests that close into 
the neutron star, where we expect high gas densities, the 
gas must have high Thomson optical depths to produce 
the blackbody-shaped spectrum of the soft excess. We 
discuss this regime in the next section. 

7. REPROCESSING BY OPTICALLY THICK MATERIAL 

Reprocessing might occur in gas that is completely op- 
tically thick, with tt > 5. Such material is naturally 
expected to reside in the accretion disk, for example, 
and when illuminated with hard X-rays it re-radiates in 
a blackbody-like spectrum. If the surface area of this 
region is of the appropriate size, it can reproduce the 
observed temperatures and luminosities. The exact ge- 
ometry of the accreting material is unknown, but for sim- 
plicity we use the phrase "reprocessing by the inner disk" 
to refer to all shapes (shells, warped surfaces, etc.) of op- 
tically thick gas at the disk's inner edge. 

To begin we consider a thick partial spherical shell 




Accretion disk 



Reprocessing 
shell 

(solid angle O) 



Fig. 9. — Schematic of the simple picture used to calculate the 
blackbody emission radius Rbb ■ 



centered on the neutron star, subtending a solid angle 
£1 at the X-ray source (see Fig. 9). Any photons that 
strike this gas are absorbed and reprocessed, so that for 
isotropic emission, 

L so ft = filoc- 

If the radius of this shell is Rbb , the reprocessing region 
has an area 

A = fi47ri?| B . 
For isotropic blackbody radiation, 

L so h = f247ri?! B crTg B . 

Thus we have 



R 



BB 



This differs from some estimates for Rbb i which assume 
all the soft luminosity is radiated from a full sphere at 
-Rbb • Our formula is more appropriate if the blackbody 
comes from reprocessing by a partially-covering, disk- 
like structure. For typical Lx — 2 x 10 38 ergs s _1 and 
kTsB = 0.17 keV, we obtain i?bb = 10 s cm. This value 
is close to the typical radius of the magnetosphere for 
the luminous XBPs, suggesting that the reprocessing gas 
may be at the inner edge of the accretion disk. 

7.1. Her X-l and Cen X-3 

There is evidence for the above picture in Her X-l and 
Cen X-3. In Her X-l, the 35-day periodicity is attributed 
to the blockin g of the central neutron st ar by a process- 
ing disk fe.g.. iGerend fc Bovntonlll976j) . Detailed geo- 
metric models for the disk ha ve been used to ex plain 
phenomena such as X -ray dips llShakura et alJl!999|) and 
varying pulse profiles ijScott et~aTl2000() . The anomalous 
low states observed in 1983, 1993, and 1999 in which the 
source remained in the low state for several 35-day pe- 
riods, have also convin cingly been shown to be due to 
occul tation by the disk ijCoburn et alJl200"ol iManchandal 
1200^ . 

Observations suggest that this blocking disk repro- 
cesses the hard X-rays to produce the soft excess. 
lEndo et all l)2000l) performed pulse-phase resolved spec- 
troscopy on ASCA data taken during the main-on state 
(</>35 ~ 0.1). As described in § 2, they found that the soft 
blackbody component and line emission features pulsed 
230° out of phase fr om the narrow power- law pulses, an 
offset seen earlier by iMcCrav et al.l l)1982j) . The authors 
proposed that the tilted inner accretion disk absorbs the 
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hard X-rays as the pulsar points away from the line 
of sight, giving rise to the blackbody and lines, which 
are emitted back toward the observer out of phase with 
the hard pulses. For the spectral parameters found by 
lEndo et alJ l|200(il) . R bb ~ 5 x 10 7 cm. They showed 
that the cooling timescale of the irradiated gas is a small 
fraction of a second, so it is possible for the reprocessed 
s oft X-rays to pulse al ong with the hard component. 

tRamsav et alJ 1)2002(1 analyzed observations with XMM 
taken at ^35 = 0.17 (close to main-on), 0.26 (low state), 
and 0.60 (short-on state). They examined how the phase 
shift between the soft excess and hard component varies 
with . They found that the profiles have phase offsets 
that vary from 150° at 35 = 0.17 to -10° at </> 35 = 0.26 
and 90° at ^35 = 0.6. They argued that this continuous 
change in offset is expected from reprocessing by a tilted 
precessing disk. 

For Her X-l, the range of observed values £;Tbb ~ 0.1- 
0.16 keV and L x ~ 2 x 10 37 ergs s" 1 gives -Rbb ~ 0.5- 
1.2 x 10 s cm. An estimate of the magnetospheric radius 
ijFrank et al.l l2002) is 

R m ~ 0.5R A ~ 1.5 x lQ % m\ n Rf' 7 B% 7 cm, 

where Ra is the standard Alfven radius. Here mi is the 
mass of the neutron star in M , Rq is its radius in 10 6 
cm, L37 is the X-ray luminosity in 10 37 ergs s _1 , and 
B12 is the neutron star surface magnetic field in 10 12 
G. For Her X-l (i? fi ~ 1 and B 12 ~4),i? m ~3x 10 8 
cm. iScott et all ll2000l) calculate R m using several other 
mo dels (i .e. lFinger et aLlll996t iKiralv fc MeszarosllT98& 
1X11119801) and find they all give R M ~ (3-4) x 10 8 cm. 
Thus, R m is larger than -Rbb by as much as a factor 
of 4 or more. We discuss possible explanations for this 
discrepancy in the next subsection. 

In Cen X-3, the soft excess is less well studied 
but is again c onsist ent with reprocessing by the disk. 
iBurderi et alJ (2000) found a bright soft excess in the 
BeppoSAX spectrum, with /cTbb = 0.11 keV and con- 
taining 58% of the total unabsorbed flux. With Lx ~ 
10 38 ergs s~\ we have R B b ~ 2.4x 10 8 cm, while B 12 ~ 3 
gives Rm ~ 1-5 x 10 8 cm. As in Her X-l, the blackbody 
and iron line are likely produced by rep rocessing; at the 
magnetosphere. iKohmura et alJ l|2001|) tested this pic- 
ture by examining aperiodic oscillations of the X-ray flux 
and measuring the time delay of the iron lines with re- 
spect to the other bands. They found a systematic time 
delay of 0.39 ±1.0 ms, which implies a light-travel dis- 
tance of ~1.7 x 10 8 cm. This distance is comparable to 
Rbb and i? m and so supports the interpretation of disk 
reprocessing. 

7.2. SMC X-l and LMC X-4 

Although SMC X-l and LMC X-4 are more luminous 
than Her X-l, they all show remarkable similarities in- 
cluding a sinusoidal soft pulse profile, and a long-term 
periodicity, which is also attributed to a blocking, pre- 
cessing disk. These similarities indicate that optically 
thick reprocessing may be the source of the soft excess 
in these systems as well. 

Given the tempera tures and luminos ities for SMC X-l 
and LMC X-4 from IPaul et alJ (|2002t see Table 1) we 
find Rbb = 1.3 and 1.1 x 10 8 cm, respectively. For LMC 
X-4 , B 12 ~ 10 indicates R rn ~ 2.9 x 10 8 cm. This gives 



Rbb/Rm — 0.4, compared with Rbb/Rtu ~ 0.2-0.4 for 
Her X-l, depending on the value for Tbb, and ~1.5 for 
Cen X-3. The values for Rbb and R m are shown for 
comparison in Table 5. 

In SMC X-l, no cyclotron feature has been detected, 
so there is no magnetic field estimate. For SMC X-l 
to have Rbb / Rm = 0.4, as in Her X-l and LMC X-4, 
would require B12 ~ 20. Such a high magnetic field 
would give a cyclotron line at ^200 keV and so may be 
difficult to detec t with current instruments. However, 
iMakishima et all |l999) suggest that the relatively low 
spectral cutoff energy (6 keV) in SMC X-l may indicate 
a smaller magnetic field than Her X-l, so the question 
is still open. Another estimate for the inner edge of the 
accretion disk is the corotation radius 

R COI = (GAZ/^g) 1 / 3 , 
which is the radius at which the Keplerian angular ve- 
locity of the disk is equal to ^ns, the angular velocity of 
the neutron star. For SMC X-l, R COI = 1.3 x 10 s cm, 
w hich is equa l to -Rbb- 

IPaul et all (2002) argue that because of discrepancies 
between -Rbb and R m , disk reprocessing may not be at 
work in the luminous XBPs. However, Her X-l and Cen 
X-3 show substantial evidence for reprocessing by the 
disk, yet have markedly different ratios between -Rbb and 
R m . It is therefore difficult to base conclusions on the 
exact values obtained for -Rbb and R rn . This uncertainty 
exists fo r several reasons, some of which are treated in 
detail bv lScott et all {20001. First, there may be errors in 
our estimate of R m , largely because the standard mod- 
els for the magnetospheric radius are incomplete. The 
interaction of the accretion disk with the magnetic field 
is poorly understood, but certainly a high rate of mass 
flow from the disk onto th e magnetic field lines has some 
effect on the field (e .g., lAgapitou fc Papaloizoul l2000t 
lAlv fc Kuiioerslll990D . Whether this makes R m larger 
or smaller is uncertain. It may also be that the inner 
edge of the disk does not actually reside at the radius 
that we define as the magnetosphere, but that some gas 
exists inside this radius. Again, such a process is not 
fully understood and warrants more investigation. 

Even if the standard formulae are accurate and the 
accretion disk really is truncated at R m , there are un- 
certainties in the parameters that are used to calculate 
R m . B is taken from measurements of a cyclotron fea- 
ture that may arise not in the main dipole field but in 
a region of stronger or weaker field near the accretion 
flow. Since R m cx -B 4 / 7 , a 25% error in the B field 
strength would cause a ~15% error in R m . There is 
also uncertainty in Lx due to errors in the distance to 
the source, uncertainty in the spectral parameters, and 
variation of the source with time. S 
25% error in Lx would cause a ^7% error in -Rm- There 
are also uncertainties in the masses and radii of neu- 
tron stars. Masses of so litary pulsars seem to be tightl y 
peaked around 1.4 M Q ijThorsett fc Chakrabartvlll999jK 
but for accreting pulsars the masses may be greater. For 
exa mple, from radial velocit ies of the companion to Vela 
X-l. lQuaintrell et all |200l obtain M = 2.27±0.17 M©. 
Since R m cx M 1 / 7 , such a mass variation would corre- 
spond to an increase in R m of ~10%. Estimates of -Rns 
vary from 0.7 to 1.6 xlO 6 cm, which can vary R m by 
another ^20%. The overall effect of these uncertainties 
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TABLE 5 

Blackbody and magnetospheric radii for selected XPBs 





Her X-l 


Cen X-3 


LMC X-4 


Vela X-l 


L x (10 38 ergs s" 1 ) 


0.2 


1.0 


1.2 


0.01 


B (10 12 G) 


4 


3.5 


10 


2.4 


kT BB (keV) 


0.16, 0.1 


0.11 


0.17 


0.2 


Rbb (10 8 cm) 


0.49, 1.2 


2.3 


1.1 


0.07 


R m (10 8 cm) 


2.9 


1.5 


2.9 


4.5 


RBB/Rm 


0.2, 0.4 


1.5 


0.4 


0.01 



is to create substantial uncertainty in R m which may be 
as high as ~30% or greater. 

There may also be errors in our estimate of Rbb- The 
spectrum of the reprocessing material is almost certainly 
not a single-temperature blackbody. A hot disk coronona 
could Compton-scatter some of the blackbody photons, 
giving an observed T higher than the true blackbody 
temperature and causing an underestimate of Rbb ■ On 
the other hand, the reprocessing surface may be tilted 
with respect to the hard X-ray flux, which would increase 
the emitting area at a given radius and thus decrease i?bt>- 
Viewing angle effects due to non-isotropic emission could 
also be important. 

Given these uncertainties, the correspondence between 
R m and Rbb to a factor of 2 or so suggests that repro- 
cessing by the inner disk can explain the spectra of LMC 
X-4 and SMC X-l. Such a model arises naturally from 
the geometry of these sources, and, unlike the other pro- 
cesses we have considered, it reproduces the pulsations 
and high luminosities of the soft component. 

7.3. RX 30059. 2-7138 and XTE J0111. 2-1311 

Transien t Be/XBP systems are expe cted to form accre- 
tion disks l|Havasaki fc Okaz aki 2004]), so inner disk re- 
processing might occur. In RX J0059.2-7138 and XTE 
J0111. 2-7317, there is no magnetic held measurement, 
so estimates of Rm are difficult. A typical value of B 
for XBPs (Bi2 = 4) gives R m ~ 1.5 x 10 8 cm for XTE 
J0111. 2-7317. This is of order the Rbb ~ 0-8 x 10 8 cm 
obtain ed from the ASCA spectrum bv lYokoaaw a et all 
(2000). RX J0059. 2-7138 has a similar luminosity and 
spectrum, and so may also have disk reprocessing. As in 
SMC X-l and LMC X-4, this is the only mechanism that 
can produce the spectral shape and high L so ft in these 
systems. 

A puzzle is the lack of pulsations in the soft component 
of RX J0059. 2-7138. One explanation is that, unlike Her 
X-l, SMC X-l, and LMC X-4, the disk may be at low 
inclination. In this case, the whole reprocessing region 
would present roughly the same angle to the observer, so 
we would see a brightly emitting disk at all phases in the 
pulsation. Another possibility is that electron scattering 
from a large disk wind or corona, of size ~10 10 cm, could 
wash out the pulsations in the soft component. 

7.4. Vela X-l and 4U 1626-61 

In Vela X-l, reprocessing by the disk is not the main 
source of the soft excess. Given Lx ~ 10 36 ergs s _1 
and B12 ~ 2.4, we obtain R m ~ 5 x 10 8 cm. Although 
a blackbody fit has not been pe rformed to the Ve la X- 
1 soft component, we note that iPaul et all l)2002|) were 



able to fit both blackbody and bremsstrahlung models 
to SMC X-l and LMC X-4 spectra, giving kT hh = 0.17 
and 0.18 keV, and kT t b = 0.3 and 0.5 keV, respectively. 
For Vela X-l, kT t b — 0.5 keV, so we assume an equiv- 
alent kT B B ~ 0.2 keV. This gives Rbb ~ 7 x 10 6 cm, 
or Rbb/Rth ~ 0.01, which is extremely small. In other 
words, Vela X-l is not luminous enough to produce a re- 
processed blackbody soft component at £;Tbb ~ 0.2 keV 
and over an area ~ i? 2 ^. On the other hand, because 
of the relatively low luminosity, the soft component is 
well explained by emission from diffuse gas excited by 
collisions and photoionization, as discussed in §§ 5-6, re- 
spectively. 

The situation in 4U 1626-67 is difficult to evaluate be- 
cause the distance to the source, and thus Lx, is poorly 
constrained. The different possibilities for a Roche lobe- 
filling companion to 4U 1626-67 imply different dis- 
tances: 1 kpc for a 0.02 Mq white dwarf, 3 kpc for a 0.08 
M Q partially degenerate, hydrogen-deple ted star, and 
36 kp c for a 0.6 helium burning star ijChakrabartvl 
Il998|) . IChakrabartvl l)1998j) found that optical photom- 
etry is consistent with emission coming mostly from an 
X-ray heated disk, which implies a distance in the range 
o f 5 kpc < D < 13 kpc. 

Orlandini ct al. (1998a) find that for a blackbody soft 
excess, the radius of the emitting region is 1.5 x D^ pc 
km. For D ~ 5 kpc and thus Lx ~ 2 x 10 36 ergs s _1 , this 
region is the size of the neutron star, so soft excess may be 
thermal emission from the star's surface. To evaluate the 
possibility of reprocessing by the disk, we take B12 = 3.5, 
kT B B = 0.3 keV, and L x = 7.7 x 10 34 L>£ pc and find that 
Rbb is of order R m ~ 9 x 10 7 cm only at D ~ 100 kpc. 
This is an unphysically large distance. We conclude that 
disk reprocessing is not at work in this system. 

7.5. Pulsations from disk reprocessing 

To conclude this section, we show that soft pulses can 
arise from disk reprocessi ng in the lum i nous XBPs. Fol- 
lowing the calculation of lEndo et al.l l)200fl) . the total 
thermal energy of the soft emitting region is simply the 
volume of the region times its thermal energy density 

Ebb = fi47r-Ri B d f ^nfcT b b^ , 

where d is the penetration depth of hard X-rays into the 
heated region. We assume that most X-rays penetrate no 
more than a few Compton depths, so that nd ~ 3/ctt ~ 
10 25 cm- 2 . For typical Rbb = 10 8 , ^/4tt - 0.1, and 
kT B B = 0.17 K, we have Ebb ~ 10 32 erg. The cooling 
timescale is simply i coo i = E B B/L so ft, which for L soft ~ 
10 37 ergs s _1 gives i coo i ~ I0 -5 s. As in Her X-l, this 
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is much shorter than the pulsation period, so soft pulses 
should be possible for this reprocessing mechanism. 

8. DISCUSSION 

We have found that only a few of the models that have 
been proposed can successfully account for the observed 
soft excess in XBPs. In addition, the origin of the soft 
excess depends strongly on the luminosity of the source. 
For systems with high L x > 10 38 ergs s" 1 (SMC X-l, 
LMC X-4, RX J0059. 2-7138, XTE J0111.2-7317, and 
Cen X-3), the soft excess can be explained only by re- 
processing of the hard X-rays by optically thick accret- 
ing material. However, for systems with Lx < 10 36 ergs 
s _1 , the soft excess is emitted by diffuse gas through colli- 
sional heating or photoionization (Vela X-l, RX J0101.3- 
7211, and AX J0103-722), or is possibly thermal emis- 
sion from the neutron star's surface (4U 1626-67, X Per). 

The more luminous XBPs probably also emit soft X- 
rays from a diffuse cloud and from the surface of the 
neutron star, but these components are dominated by 
the flux from the power law and from disk reprocessing. 
On the other hand, in lower Lx systems we do not expect 
disk reprocessing to be important, since for blackbody 

emission at Rbb = Rm, we have Tbb oc L^/ 4 / Rm 2 . Since 

R m oc L x 2 ^ 7 , this gives Tbb oc L^ 28 . Thus for lower 
Lx the temperature of the disk moves toward the EUV 
and out of the soft X-ray band. 

For XBPs with moderate Lx ~ 10 37 ergs s _1 , it ap- 
pears possible to have disk reprocessing (Her X-l), emis- 
sion by diffuse gas (EXO 053109-6609.2), or possibly 
both (A 0538-66). We thus have something of a con- 
tinuum, in which the soft excess from more luminous 
XBPs comes primarily from disk reprocessing, while in 
less luminous ones it originates in optically thin emission 
or thermal emission from the neutron star's surface. 

In addition to these general conclusions, this work 
raises a number of puzzling questions. One is the appar- 
ently high brightness temperature of the power-law com- 
ponent in SMC X-l (§ 4.1). This may suggest that the 
accretion column is hotter or larger in extent than mod- 
els predict. For Vela X-l, the question remains whether 
the soft excess emission is powered mainly by photoion- 
ization or collisional energy sources (§§ 5-6). In RX 
J0059. 2-7138, where disk reprocessing is likely at work, 
the complete absence of soft pulsations (§ 7.3) is also 
worth further investigation. 

It is important to examine the lack of exact correlation 
between Rbb and R m in the sources where we believe 
reprocessing by the accretion disk is at work. We have 
discussed a number of uncertainties in these calculations 
in § 7.2, but this problem warrants further study. For 
example, a thorough investigation of non-blackbody ef- 
fects in the disk reprocessing can help reduce uncertainty 
in Rbb ■ Radiative transfer in an obliquely il luminated 
model disk has been treated bv iPsaltia l)2002j) . and the 
structure and line emission of an i lluminated disk atmo- 
sphere has been modeled in detail l)Jimenez-Garate et alJ 
120011 120021) . Using such models to analyze the soft emis- 
sion from XBPs can help diagnose the gas conditions near 
the disk's inner edge. 

There remains the question of how the accreting ma- 
terial, which is expected to be in a geometrically thin 
disk, is able to subtend a substantial solid angle to the 



central neutron star. Interactions with the magnetic 
field are likely t o be key in this process. For example, 
iSpruit fc Taaml l)199QD have shown analytically that ac- 
creting matter can be tied to the magnetic field in equi- 
librium positions off the midplane of the disk, near its in- 
ner edge. Numerical, three-dimensional magnetohydro- 
dynamic simu lations of accretion onto rotating magne- 
tized stars llRomanova et al.ll2003l 1200*2]) show structures 
including funnel flows and dense "rings" that develop 
near R m . These may be the sites of the X-ray reprocess- 
ing, and the soft excess emission may provide constraints 
on the properties of such structures. 

We can also use the soft excess to study warp- 
ing and precession of the disk, particularly in Her X- 
1, SMC X-l, and LMC X-4. A warping instability 
in a centr ally-illuminated disk was found theoretically 
by iPringld (l996j) a nd has been exp l ored i n a num- 
ber of studies ("e.g.. jWiiers & Pringlel 119991 IL"al Il999t 
iMalonev fc Bege lman||19 97fl that also predict the disk's 
precession. As iRamsav et all {2002) have done in the 
case of Her X-l, it would be interesting to examine SMC 
X-l and LMC X-4 and see how the reprocessing region 
on the disk moves as the disk precesses. 

9. SUMMARY 

In this paper we have reviewed a number of observa- 
tions of the soft excess in X-ray pulsarsand explored the 
origin of this feature. We find the following: 

1 . The soft excess is a very common, if not ubiquitous 
feature of emission from XBPs. Observations of 
XBPs that have not shown a soft excess are most 
often limited by low flux, high absorbing column, 
or insufficient soft sensitivity. 

2. For luminous (Lx > 10 38 ergs s _1 ) sources, the 
soft components have high luminosities and soft 
spectral shapes that can only be explained by re- 
processing of hard X-rays from the neutron star by 
optically thick, accreting material, most likely near 
the inner edge of the accretion disk. 

3. For less luminous (Lx < 10 36 ergs s™ 1 ) sources, 
the soft excess is due to other processes, such as 
emission by photoionized or collisionally heated dif- 
fuse gas, or thermal emission from the surface of 
the neutron star. In these sources the soft compo- 
nent cannot come from reprocessing by the accre- 
tion disk. 

4. XBPs of intermediate luminosity (Lx ~ 10 37 ergs 
s _1 ) show either or both of these types of emission. 

5. Future observations of soft excess will allow us to 
explore aspects of XBPs such as the properties of 
circumstellar gas, interaction of accreting matter 
and the magnetic field, and warping and precession 
of the accretion disk. 
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